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Abstract 

We study the 0{N) symmetric linear sigma model at finite temperature as the 
low-energy effective models of quantum chromodynamics(QCD) using the Cornwall- 
Jackiw-Tomboulis(CJT) effective action for composite operators. It has so far been 
claimed that the Nambu-Goldstone theorem is not satisfied at finite temperature in 
this framework unless the large limit in the 0{N) symmetry is taken. We show 
that this is not the case. The pion is always massless below the critical temperature, if 
one determines the propagator within the form such that the symmetry of the system 
is conserved, and defines the pion mass as the curvature of the effective potential. 
We use a new renormalization prescription for the CJT effective potential in the 
Hartree-Fock approximation. A numerical study of the Schwinger-Dyson equation 
and the gap equation is carried out including the thermal and quantum loops. We 
point out a problem in the derivation of the sigma meson mass without quantum 
correction at finite temperature. A problem about the order of the phase transition 
in this approach is also discussed. 
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1 Introduction 



Chiral symmetry is one of the most important features of low-lying liadron properties. In 
quantum chromodynamics(QCD) this symmetry is well satisfied in the SU (2) sector due to 
the light u, d quark masses. Existence of the nearly massless mesons, or pions, means that 
the chiral symmetry must be spontaneously broken and the Nambu-Goldstone(NG) bosons 
appear. Spontaneous chiral symmetry breaking(SCSB) is also manifested in the non- 
degenerate parity doublet of nucleons. So SCSB influences the low-lying hadron spectra 
significantly. Though such low-energy phenomena of the strong interaction as well as the 
quark confinement are now widely confirmed in experiments, no analytic investigation 
based on the first principle, or QCD is carried out due to its non-perturbative nature. So 
it is still important to study the effective models of QCD in order to understand the nature 
of hadrons in addition to numerical analyses in the lattice gauge theory. 

In increasing temperature, as shown in the lattice QCD calculation, the chiral symmetry 
is believed to be restored at around T = 100 ~ 300 MeV. The study of physics at finite 
temperature is very interesting from both theoretical and experimental points of view. 
The big bang model tells us that a series of phase transitions, which of course include the 
QCD phase transition, occurred in the early universe. It could also be possible to probe 
the underlying physics of QCD in laboratory involving relativistic heavy ion collisions. 
These experiments are planned in near future and its result is expected to elucidate some 
important questions like the mechanism of the chiral symmetry restoration and the nature 
of quark gluon plasma. 

Theoretical investigation of the symmetry restoration at finite temperature in terms of 
field theory was first studied by Kirzhnits and Linde[0. They observed that spontaneous 
symmetry breaking(SSB) will be restored at sufficiently high temperatures. It is now well- 
known that the naive effective potential up to the 1-loop level does not work at finite 
temperature. Weinberg pointed out that at very high temperature, powers of temperature 
T can compensate for powers of a small coupling constant e, leading to a breakdown of the 
perturbation expansion|0. He showed that the leading effect of this sort arises from the 
e^T^ term in the scalar 0"^ theory. Dolan and Jackiw and others showed that systematic 
summation of a certain kind of loop diagrams in this model is needed at finite temperature 
ii- 

For analysis of the linear sigma model at finite temperature, the mean field approx- 
imation has mainly been used so far. The mean field theory is a theory in which the 
resummation of the diagrams is introduced approximately by hand and has been used by 
many authors ||, 0] . 

Recently several analyses at finite temperature have been done using an extended effec- 
tive potential for composite operators introduced by Cornwall, Jackiw, and Tomboulis(CJT) 
Contrary to the usual effective action, their effective action depends not only the 
classical field (pc^x) but also on G{x,y). These two quantities are to be realized as the 
expectation values of a quantum field and the time ordered product of the field op- 
erator T(f){x)(j){y) respectively. In this case the effective action r(0c,G) is the generating 
functional of the two particle irreducible vacuum graphs. This formalism was originally 
written at zero temperature but it has been extended for finite temperature in cf)^ theory 
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by Amelino-Camelia and Pi0. 

There is an advantage in using the CJT formahsm to calculate the effective potential in 
the Hartree-Fock approximation. According to ref.0, we need to evaluate only one graph 
that of the "double bubble" instead of summing infinite many "daisy" and "super daisy" 
graphs using the usual tree level propagators. These infinite diagrams are incorporated 
automatically in the effective action if the CJT formalism is used, so the effective potential 
takes real values at finite temperature. The need of the resummation of this kind of loop 



diagrams at finite temperature is also discussed in other approaches |]T0 . 

Amelino-Camelia et.al. investigated the 0{N) linear sigma model, which is regarded 
as the low-energy effective model of the sigma meson and the pions, at finite temperature 
with the CJT action^, It is known, however, that in their approach one encounters 
a difficulty that the NG theorem appears to be violated in the broken symmetry phase at 
finite temperature. They concluded that in the CJT formalism the NG theorem is violated 
in the Hartree-Fock approximation, i.e., the case of finite values of and the so-called 
super-daisy approximation, while it is satisfied if and only if leading contributions in the 
1/A^ expansion are taken. They also discussed the renormalization in the CJT formalism 
and suggested that a consistent renormalization cannot be performed in the broken phase 
since the non-perturbative quantities, or solutions of the Schwinger-Dyson(SD) equations, 
are included. 

In this paper we study two subjects. One is the formulation of the linear sigma model 
at finite temperature. We re-analyze the 0{N) linear sigma model at finite temperature 
in the CJT formalism and show that the NG theorem is always satisfied at any finite 
values of A^ in the Hartree-Fock approximation. The CJT action does not violate the 
0{N) symmetry, so that the pions are always massless below the critical temperature 
if we define the meson masses by the curvature, or the second derivative of the effective 
potential. On the renormalization, we adopt the method of the renormalization of auxiliary 
fields. As shown in the CJT's paper, the solutions of the SD-equation can be thought as 
a kind of auxiliary field in the mean field level. We show that the CJT effective potential 
is renormalizable if the technique of the auxiliary field is applied. This approach has a 
clear advantage that once the effective potential is renormalized, all the physical quantities 
derived from the effective potential are finite. 

The other is the application of the above formalism to low-energy hadron properties. 
The linear sigma model has various merits as an effective model of low-energy hadron 
dynamics. It can describe SSB as in QCD. Furthermore one can study both the symmetry 
broken phase and the restored phase if desired. Other effective models such as non-linear 
sigma model or chiral perturbation theory treat basically the world with SSB. 

The contents of this paper are as follows. In Sec. |^, we introduce the CJT effective 
action for composite operators and explain the difference from the ordinary effective action. 
In Sec. ^ we construct an effective potential for the 0{N) symmetric linear sigma model. 
Here we also point out defects of the conventional formulation. We show that the NG 
theorem holds even at finite temperature if one adopts the 0{N) symmetric form of the 
propagator and defines the meson masses as the curvatures of the effective potential. We 
also give a new renormalization prescription of the CJT effective potential. It is based on 
the renormalization of auxiliary fields and was long ago applied to this model in the large 
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limit [|1^, [T^. It can remove divergences in the effective potential both in the symmetry 
broken and unbroken phases, while the conventional renormalization can be performed only 
in the restored phase. Using these techniques, we apply the linear sigma model to the low- 
energy mesons in Sec. |[ As is well-known, the 0(4) linear sigma model can describe the 
physics of the sigma meson and the pions in low-energy region. We discuss temperature 
dependence of physical masses of these mesons and sigma condensate corresponding to 
the pion decay constant. In Sec. we concentrate on the order of the phase transition 
and consider the so-called setting-sun diagram. We discuss approximate treatments of the 
setting-sun diagram and compare them. Finally summary and conclusion are given in Sec. 

0- 



2 Effective Action for Composite Operators 

Effective action for composite operators was originally studied in condensed matter physics 
p!6| and then extended to relativistic field theories . We here use an approach introduced 
by Cornwall, Jackiw and Tomboulis [^. This is based on functional methods or the path 
integral representation for the Green functions. 

The effective action for composite operators is a generalization of the conventional ef- 
fective action and is written by r[0c, G]. This is a functional both of the expectations value 
of the quantum field (pdx) = {0\(j){x)\0) and of the propagator G{x,y) = (O|T0(x)0(?/)|O). 
The c-number function (pd^) is also called a classical field. The variational equations 



0, (1) 
= (2) 



6G{x,y) 

determine 0c and G on the vacuum. The equation (Q) is nothing but the SD-equation 
for the propagator G. Furthermore, we can show that the second variational derivative 
of r[0c, G] in G leads to the Bethe-Salpeter equation which describes relativistic bound 
states. For example, in hadron physics, this is used in order to describe meson states as 
the bound state of the quark-antiquark pair. 

We now describe the series expansion for r[0, G]. We introduce the so-called tree level 
2-point Green function by 

tV {(j)c,x,y) = =iD X - y + 3 

d(t)c[x)b(t)c[y) d(pc{x)d(f)c{y) 

The required series is then 

r[0„ G] = S[<f),] + ^TrLnG-i + '-TtV-\<P,)G + Fal^c, G] + const (4) 

where the trace, the logarithm and the product V~^G are taken in the functional sense. 
The conventional effective action F[0c] is F[0c, G] with / = 0, i.e., 

F[0e] = r[0c,Go] (5) 
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^0^.0 (6) 

6Go{x,y) 

The constant which is independent of and G is evaluated so that eqs.(||) and (||) are 
satisfied: 

r[0„ G] = S[^c] + '-TrLnDG-' + '-Ti(V-\(l>,)G - 1) + T^icPc, G] (7) 

This expression is used in some hterature. r2[0c,G'] is given by all the two-particle and 
higher two-particle irreducible vacuum graphs in a theory which has the vertices determined 
by the interaction of the action 5'int[0, 0c] and the propagators G{x,y). 

When one considers the case of translation- invariant solutions, one sets (pdx) to a 
constant 0c and takes G{x, y) to be a function only of x — y. The series for the effective 
potential V{(f)c, G) can be easily obtained from eq.(H): 

V^(0c, G) = Vo(0c) - ^ / ^lnG-i(fc)-^/^trI)-i(0„fc)G(fc) + V2(0c,G) (8) 
where Vo(0c) is the tree-level (classical) effective potential, 

Gik) = [ d^xe^'^'^-y^Gix - y) (9) 



l)(0c, k)= j d^xe'^-^^-y^V{(t>c. x-y) (10) 

and — V2(0c, G) is the sum all two and higher order loop two-particle irreducible vacuum 

graphs of the theory with vertices given by 5'int(0c) and propagator G{k). The stationary 
requirements are then 

^vf'^K, (u) 

<90c 

^n^^^G) = (12) 

dG{k) ^ ^'^^ 

These equations are the starting point of our discussion. 



3 Formulation of the 0{N) Linear Sigma Model 

In this section we construct the effective potential of the 0{N) symmetric linear sigma 
model using the CJT formalism. In QCD with two-flavor massless quarks, the chiral 
SU{2)l X SU{2)ji symmetry is satisfled in the lagrangian level. But this symmetry is 
spontaneously broken to the SU{2)v symmetry, because the vacuum, i.e., the ground state 
of the fleld conflguration violates it. We can see this phenomenon in the 0(4) linear sigma 
model, because the 0(4) symmetry has the same algebra as SU{2) x SU{2). Four flelds 
in the 0(4) linear sigma model are identifled with one sigma meson and three pions when 
we regard it as an effective model of hadrons. By spontaneous symmetry breaking from 
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0(4) to 0(3), the three fields become massless according to the NG theorem and the 
remaining one is massive. Here we concentrate on the meson fields and see how this model 
is formulated with the CJT effective action. We set the arbitrary N instead of = 4 at 
the stage of the formulation for generality. We use the Euclidean metric hereafter because 
it is convenient when we extend the model to that with the imaginary time formalism at 
finite temperature. 

The lagrangian density for the 0{N) linear sigma model is given by 



2^ 

6iv' 



(13) 



where 0^ = (f)°'(j)°', a runs over 1 to and repeated indices are summed. By shifting the 
field as 0"(x) = (j)"'{x) + The tree-level 2-point Green function V is obtained: 



4A 

6^'6\x-y) + —rc<l>y{x-y) (14) 



In momentum space, it is given by 
and the inverse becomes 

1 



k"^ + + —c 

3N' 



+ 



3N 



1 



(15) 



i;ab 




The interaction lagrangian which describes the vertices of the shifted theory is given by 



(16) 



2A 

3iv' 



+ 



6iv' 



(17) 



The diagrams contributing to r2 [00 G\ are shown in Fig. |l|. Each line represents the 
propagator Gab{x,y), and there are two kinds of vertices: a four-point vertex proportional 
to A and a three-point vertex, which results from shifting the fields, proportional to X(j)'^{x). 

If we were computing the ordinary effective action r[0c], the lines would represent the 
propagator T>ab{4'c, x, y) and there would be additional contributions which are two-particle 
reducible. These diagrams are shown in Fig. 0. 

Let us compute the CJT effective action(|). We evaluate r[0c-O] in the Hartree-Fock 
approximation, i.e., we take into account the "cxd" type diagram only. This diagram is a 
leading order in r2[0c;O] in both the loop expansion and the 1/A^ expansion. The other 
two-loop diagram, the setting-sun type diagram, which is the next to leading order in the 
1/A^ expansion, is quite lengthy to calculate and causes a problem on renormalization (see 
later) . 
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In this approximation, r2 is given by 

A f 

V2{^c-,G) = — J d^x{Gaa{x,x)Gbb{x,x) +2Gab{x,x)Gba{x,x)) 

Then the CJT action has the form 

r[0c,G] 



(18) 



+ ^TrLnG-i + ^Tr{V-\(j),)G) 



A 



6N 



d^x{Gaa{x, x)Gbbix, x) + 2Gah{x, x)Gba{x, x)) 



where S[4>^ is the classical action 

s[4>.] = / d'x {^-d,rA€ + \m'4>l + ^{<t>l? 

From this, the effective potential for the composite operators is given by 
V{<t>c, G) 

= + IJ -^,\ndetG~\k) + \ j ^MT^-\'i>c,k)G{k)) 



(19) 



(20) 



A 



{k) 



d^k 



Gba{k) 



with 



(22) 



where the trace and determinant are taken only in the internal (flavor) space. 

If one applies this model to a system at finite temperature, one replaces the loop integral 
with 

d^k 1 A /• rf^fc 



(2vr) 



E 



i 

T(3,k 



(23) 



in the imaginary time formalism. Here /? = 1/T is the inverse temperature and the sum is 
taken over the Matsubara frequency. So the starting effective potential we study is 



(24) 



Minimizing the effective potential with respect to the propagator G{(l)c,k), we obtain 
the SD-equation 

2A ,^ „ „ ... ^25) 



GaM = Kbi<Pc,p) + ^ [f.^G,c{k)6'^' + '^i^Ga,{k) 

The solution Gq{4>c-, k) of this equation is inserted back into the expression for the effective 
potential to give the conventional effective potential as a function of 0c- 
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3.1 How to Take the Form of the Propagator G 



3.1.1 Conventional Forms of the Propagator G 



First we need to determine the form of the SD-equation, G. In case of the 0{N) linear 
sigma model, it seems to be natural to adopt the following multi-component form, 



1 ^ 1 

a J — z TV 



(26) 



In fact several authors adopt this form||l2|, [13[ and in ref.|jT3| this is called "dressed propa- 
gator" ansatz. They also identify the parameters and with the masses of the sigma 
meson and pions, respectively. It is, however, obvious that this form of G does violate 
the 0{N) symmetry. So the effective potential after substituting the solution Go also vi- 
olates the 0{N) symmetry. This is one of the reasons why this formulation erroneously 
concludes that the pions which should be massless below the critical temperature acquire 
finite masses. 



3.1.2 An 0{N) Symmetric Form of G 

We calculate the meson masses by using the definitions (0) and (^) from the effective 
potential directly. In the CJT formalism, we need to choose the form of G not to break 
the symmetry. If one takes eq. (p6D , for instance, we will lose the 0{N) symmetry in the 
resulting effective potential. Here we use the following form of G from the analogy of V. 



ab 



+ 



Y( 



N 



(27) 



where X{(j)l) and Y{(f)'^) are unknown functions which are determined by solving the SD- 
equation. We assume that X and Y are independent of the external momentum and also 
that no more parameters are needed because the Hartree-Fock approximation is employed. 
Substituting eq.(^) into the CJT effective potential(|D, we obtain 



y(0„G(0„A;)) 



,9 / ,9x9 1 ^r, ^ 1 X 



-g[M 



x\ 



1/9 2A 

-- H < 

2 I N 



Mly P[Mxy] + 



A^- 1 



m + 



2^ ^2 



Mi 



+ 



A 



where we introduce 



W[MxYf + (A^^ - l)P[Mxf + 2(A^ - 1)P[Mxy]P[M_ 



x\ 



X 



X( 



N ^ 
N 



P[M_ 



x\ 



(28) 

(29) 
(30) 
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and 



Q[m] = ^ ln(A;^ + m^) 

Tf3,k 



(31) 
(32) 

°(/>" and therefore is 0(-/V) 



for simplicity. This effective potential depends only on (pl i 
symmetric obviously. 

The functions Mx and Mxy are determined from the following SD-equations. 



M'xY 



Mi 



m + 



2A 



1 



2A 



AT 



2AA^+ 1 



P[M 



AT 



(33) 
(34) 



Using these equations, the gap equation for the sigma condensate is given by 



dV{4>c,G{<Pc)) 



dan 



( dV{ct>,,G{(t>c)) ^ dV G{(t>c))dG{ct>^ 



<Tc=<yo,'I^c=^ 



da, 



dG 



dar 



crc=o-Q,-rTc=0 



aV(0e,G(</.e)) 



O'c=cro,7rc=0 



which gives 



2A 



2A 



^ + ^'^o + -^P[Mxy] + 



2XN-1 







This equation reduces to the very simple form when the SD-equations are used again 

4A 



Mly- 



2,N 







(35) 
(36) 

(37) 



Here we need some case for the meson masses. The question is whether Mo- and Mn-, 
which are the parameters that appear in the SD-cquation, can be identified as the physical 
masses of the sigma and pions. In fact, it is more appropriate to define the masses by the 
curvatures of the effective potential around its minimum, namely. 



2 <^^^ 



('■c=O'0,Tc=0 



d^V 



d-Kl 



(38) 



(39) 



rather than and M^. The quantities M^. and are merely regarded as variational 
parameters which are determined from the SD-equation. In fact if higher order loop cor- 
rections are taken into account in r2 in the CJT action, the solutions of the SD-equation, 
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Mo- and M,r, become momentum-dependent and then cannot be simply identified with 
the physical meson masses. Of course, the masses defined by eqs. (pS]) and (|39|) are also 
approximate values unless all the loop contributions are taken into account. But we see 
that the NG theorem is always satisfied because the 0{N) symmetry is satisfied at any 
order of expansion. This way of defining the masses in the effective potential formalism is 
also seen in ref . 0, |18| 0. The differences between two definitions come from the truncation 



of the loop expansion of the CJT action. Since the physical mass is defined by the pole 
of the two-point function, if one can calculate the effective potential up to all-order, both 
definitions might give the same mass. Our definition of the mass is based on the approxi- 
mation, which respects the 0{N) symmetry and therefore the NG-theorem is satisfied at 
each loop diagram of r2 in the CJT action. 



We evaluate the meson masses according to the definitions (pSQ and (|39|) . The results 

are 

4A 



,2 o,2 I dMxY 



m„ = 2cr, 



rl(^2 



dm 



n SAT 



(40) 



= (41) 

in the broken symmetry phase. The pion mass vanishes as expected. The derivative term 
in the sigma meson mass is transformed into 



dMly 



X (qA^^ + 6AA^(A^ + l)P2[Mx] + 18XNP2[Mxy] 

+8{N + 2)X^P2[Mx]P2[Mxy]Y' (42) 



with 

by the SD-equations. 

Finally we comment on the differences between the conventional method and ours again. 
As we showed previously, the CJT effective potential has two arguments, i.e., 0c and G. 
Then the chiral symmetry breaking may appear in two ways, one as a nonzero 0c and 
the other as a nonsymmetric G. The first one corresponds to the ordinary spontaneously 
symmetry breaking and to our parametrization of G. The latter case is that the vacuum 
expectation value of the two point function breaks the chiral symmetry, and the dressed 
propagator ansatz(PBD corresponds to this choice. Since there occurs spontaneous sym- 
metry breaking in both cases, the NG bosons must appear. In the latter case, however, 
they do not correspond to the pions, as they are unphysically massive. If one is interested 
in low-energy QCD phenomenology or SCSB, the former choice is the right way of the 



"^In Q , the dressed propagator ansatz (g6|) are used. In this case even if one defines the meson masses 
by the second derivatives of the effective potential, the NG theorem is violated in the CJT formalism 
because of the violation of the 0{N) symmetry in the effective potential. 
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symmetry breaking because the pions appear as the massless NG bosons, and the ordinary 
effective potential is reproduced after ehminating G. The difference between the two ap- 
proaches might come from the truncation of the loop expansion of the effective potential, 
for both the approaches should give the same physical mass in the full-loop calculation. 

We discuss the treatment of these equations in a phenomenological application in Sec. 
^ together with numerical calculation. 



4 Renormalization 



4.1 A Conventional Cut-ofF Scheme 



Renormalization of the CJT effective action in the linear sigma model is discussed in 
ref.[|T5| in the past. In that paper, divergent integrals are regularized with a cut-off and 
the renormalized mass and coupling are defined in order to absorb the divergent part. This 
renormalization program was performed in the large N limit ||14|| , where all the divergent 
terms are absorbed into the bare quantities. However, if one applies this to the present 
case, one sees that not all the divergences can be removed away. We here review this 
renormalization following Amelino-Camelia [O] and point out what is the problem. 
We write the effective potential and the SD-equations for eq. (|2! 




2A 



m 



2X 2 



3N 
2A 

~ 3iV 



+ Iq[m.] 



N -1 



2A 



1 



2N 



6 N 



..2 2 2A 2 2A , 2AA^- 



N 
P[M^] 



Ml 



2 2A 2 2A + 1 , 2A , 

3A^ ' 3 ^ ^ 3A^ ^ ^ 



(44) 

(45) 
(46) 



where we take the conventional 0{N) violated form of G according to |jT2|, but our discus- 
sion on the renormalization is unchanged if we take the 0{N) symmetric form of G. These 
quantities are all written by bare quantities. The loop integral P[M], shown in eq. (PT[) , of 
course diverges. We regularize the temperature independent term with the cut-off A. 



P[M] = h- M^l2 + Pf[M] 



PfiM] 



In^r- 

/i2 



167r2 



(2^ 



exp 



A;2 + M2 



T 



(47) 
(48) 



11 



h 



A' 

8^ 



2' 



167r2 



In- 



A2 



(49) 



where /i denotes the renormahzation scale and is arbitrary. Divergence also appears in the 
integral of logarithm in the effective potential. 



\q[m] 



-—l2 + —h + Qf[M] 



QAM] 



In — + Tin 



647r2 



exp 



V 



A;2 + M2 



(50) 



(51) 



After the regularization, one must renormalize them. Amelino-Camelia determines the 
renormalized mass and coupling as 



m 

X 



R 



m 

T + 3 



2N + 2 



N 



1 



1 
A 



respectively. This is, however, renormalized only at Mg- 
M is hold, we have 



2N + 2 
3 N ' 



A£r. In fact when M^- 



(52) 



m 



3 N (^^' 
However, if we use it in the broken phase, all the terms in the SD-equation and the 
effective potential are expressed not only by the renormalized quantities but bare quantities 
remain in several terms. Essentially the same situation is encountered in the mean field 
approach [0. 

Under these circumstances, Amelino-Camelia makes the following arguments. If we are 
interested in the possibility of renormahzation, we drop the A term in the infinite limit of 
the cut-off. We then get, however, M„ = M^^. He also says that if one treats the model as 
the low-energy effective model, the cut-off A should be a large finite value which is smaller 
than the Landau pole. He concludes that vanishes at A ^ but remains finite for the 
finite value of the cut-off. 

We agree that the 0{N) linear sigma model in four dimensions may be trivial if the 
full quantum effects are included. Our analysis is however basically one-loop level and 
therefore semi-classical. So we consider another way of the renormahzation which removes 
divergences order by order both in the symmetry restored phase and broken phase. 

Recently a regularization scheme, which is similar to the Pauli-Villars regularization, is 
used in the C JT action pO| . We here give a new renormahzation scheme based on auxiliary 



field method with the dimensional regularization. 



4.2 An Application of the Auxiliary Field Method 

In order to solve the above problem of renormahzation of the CJT effective action, we here 
use the renormahzation procedure of the auxiliary field. In the original CJT's paper [||, it is 
reported that the effective potential in the large N limit in the CJT formalism is the same 
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form obtained from the auxiliary field method [p^. Since the solutions of the SD-equation, 
X and Y, are momentum independent here, we can regard them as an auxiliary field. Then 
one can renormalize the CJT effective potential with the auxiliary field method following 



some literature 14, 15 



Following the above argument, we first rewrite the 0{N) symmetric CJT effective 
potential (p8D using the SD-equations(|33|) and (|34D 



^ .4 3A^ ..2 2 3iV(iV- 1) , . 



V{(f)„MxY,Mx 

.2 A 



.^^^^{(iV+l)M^, + 3(iV-l)M^ 



+- 



-2(A^ - 1)M^yMx + 2Nm^} 
d^k ^ ,,2 ,.2 ^ N-1 r d^k 



(27r 



(2 



(54) 



Throughout this section, we ignore the temperature dependent terms since they are all finite 
and not renormalized. This means that the loop integrals in the equations are simply four- 
momentum integrals. We consider that the extension to the system at finite temperature 
should be carried out for the well-defined effective potential which has no divergence after 
renormalization. Once the renormalization is carried out at T = 0, it is straightforward to 
extend it to the system at finite temperature since the temperature dependent terms all 
converge in the UV limit due to the Boltzmann factor. 

We add four counter terms in order to cancel the divergence. 



A^- 1 



AMlym^ + B{N - l)M]^m^ + -CM^y + ^ 



-DM 



X 



(55) 



where A, B, C and D are renormalization constants to adjust the divergences and become 

1/2 



A = B = 0, 



C = D 



327r2 \A-d 



7s + ln(47r; 



(56) 



if the MS scheme of dimensional regularization is used. These correspond to the following 
renormalization conditions. 



dV 



dMly 



dV 



a2=0,A4y=M2=0 



dMl 



d^V 



{dM?--^^ 



3Nm^ 
~ 4X{N + 2) 

3N{N-l)m^ 
4A(iV + 2) 

3A^(A^ + 1) 
~ 8X{N + 2) 



(57) 

(58) 
(59) 
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respectively. Therefore the renormahzed effective potential is given by 



^N{N - 1) 
8A(iV + 2) 



2 ^^^^ SAT^'^ 4A(iV + 2) 4A(A^ + 2) ^ 

O AT 

-2(A^ - l)M|-yM|. + 2Arm^} 
647r2 2; 647r2 I ^2 



ta^-S +^ 'n^-S (61) 



The renormahzed SD-equations are the same form as eqs. (|33D , (|3^ ) except for the divergent 
terms. 

,.2 2 2A . 2XMIy f, Mly A 2\N-lMl (Ml \ 

^ SA^^'^ 3A^ 167r2 \^ /i2 ^ ^ 3 A^ 167r2 \ I 

The gap equation for the sigma condensate and the formulae of the meson masses are the 



same form as eqs.(p7D, (^OD and (|41|). They are all finite because Mxy and Mx are finite. 



Here we make a short remark on the problem in the large A^ limit. In refs.[|T5], [19 1 
the stability of the effective potential of the linear sigma model in the large A^ limit was 
discussed. They found that the true minimum of the effective potential gives the trivial 
phase in the large A^ limit. Although the renormalization scheme they use is different from 
ours, there exists the same situation in the present approach. In fact one can show that 
their results are identical to ours in the large A^ limit. Also in this limit, Mxy vanishes and 
the renormalization conditions and correspond to eqs.(2.5) and (2.6) in ref. |T5 
respectively. 



4.3 Summary 

We summarize the results of this section. First we formulate the 0{N) linear sigma model 
with the CJT effective action. The Hartree-Fock approximation is used, which corresponds 
to incorporating only the double bubble diagram in the two-loop level in the CJT formalism. 
We adopt the 0{N) symmetric form of the propagator G, otherwise the NG theorem is 
violated. On the renormalization, we use the method of the auxiliary field. One can remove 
all the divergence in the effective potential even in the broken phase in this method. The 
effective potential is then eq.(|6T|) plus the thermal effects. 

V{<f),,MxY,Mx) 
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A 

3iv' 



+ 



A\{N + 2) 



A\{N + 2) 



16A(A^ + 2; 



{(iV+ l)M^y + 3(A^- 1)M; 



-2{N - 1)M]^yMx + 2Nm^} 



+ lQTmy] + ^ 



647r2 



In 



with 



2^ 



A^- 1 
647r2 



In- 



ln(l 



(27r)3/5 

Mx and Mxy are determined by the SD-equations 



2 2A ,2 2A„ , 2AiV- 1 „ 



Ml 



N 167r2 
2A 



2\N -I Ml 
1 167r2 



2A M- 



XY 



3N 167r2 



In 



with 



PtH = J 



3 167r2 



(27r) = 



(64) 
(65) 



(66) 



(67) 



(6J 



The sigma condensate is given by eq.(p7D which is finite since Mxy is finite in this case. 

The meson masses are also given by eqs.(^OD and (|41|). The pion mass is of course 
massless in the broken symmetry phase. In the equation of the sigma meson mass, P2['^] 
factor is replaced by 



P2[m] P2,T['m] 



167r2 



In 



m 



P2,tW\ 



(2^^ 



1 + 



(69) 
(70) 



in the MS scheme of the dimensional regularization. 



5 An Application to the Low-energy Mesons 

In this and the following sections, we apply the CJT formalism introduced in Sec. ^ to 
systems of low-energy hadrons. We include both the thermal and quantum corrections in 
numerical calculation. 
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5.1 Formulation 



First we formulate a system with a sigma meson and pions. Since the real pions have 
masses because of the small but non-zero u, d quark masses, we introduce a term which 
explicitly breaks chiral symmetry in the lagrangian. 



C 



+ 



2^ 

6iv' 



k'2\2 



— CCT 



(71) 



with (j)^ = a. Since the two-point function is obtained from the second- derivative of the 
lagrangian, this term does not affect it and the SD-equations for Mxy and Mx are the 
same forms as eqs.(^) and (l67|). Thus the effective potential is the same as eq . (|6^ except 
for the CO" term. 



V{(t>c,MxY.Mx] 



M'xy4>1 
3N 



2^ 



— ca 



AX{N + 2) 

3N 
16A(A^ + 2) 



Mlym^ + 



4A(iV + 2) ^^^^"^ 



+ l)Mj^y + 3(iV - l)M% 
-2{N - 1)M]^yMx + 2Nm'^] 



+ 



In 



M'xY 



+ 



N-1 
647r2 



647r2 y /i^ 

The gap equation for the sigma condensate in turn becomes 

4A 



In- 



M 



X 



3 



3A^ 



-On 



(72) 



(73) 



From this equation one sees that the pion mass does not vanish even below the critical 
point 

(74) 



c 



and the sigma meson mass is given by 



2^0 



O"c = O"0i''''c = 

where the derivative term is the same form as eq.(^ 




(75) 



5.2 Numerical Results 

Here we show some numerical results of the above equations. We first concentrate only on 
the thermal effects and ignore the quantum corrections. This approximation is frequently 
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used in literature when one is interested in the thermal effects. We will include the quantum 
corrections in the next subsection. Throughout this and the next section, we set = 4, 
i.e., the SU{2)l x SU{2)ji chiral symmetry which corresponds to the real world of one 
sigma and three pions. 



5.2.1 Thermal Corrections 



We use the values at zero temperature as initial conditions of numerical calculation. Be- 
cause we ignore the quantum corrections at this stage, the solutions of the SD-equations 
MxY and Mx are equal to mo- and m^r at T = as in the chiral limit. Initial parameters 
are determined from the condition at T = 0: 



c=UmliT = 0) 



(76) 



A = -K(T = 0)-m^(T = 0)) 



m 



J TT 

--ml(T 



0) + ^ml{T 



0) 



(77) 



(7J 



where /,r(= 93 MeV) is the pion decay constant at zero temperature. 



For the mass of cr. Particle Data Group gives the value of 400-1200 MeV[21]. We take 



600 MeV as typical values. Qualitative properties change little if the mass of the 



sigma meson is to be rrio 
breaking takes place. 



1 GeV. m is negative because the spontaneous symmetry 



These results are essentially the same as ref . [|T3| because we ignore the quantum correc- 
however, the author identifies Mxy and Mx with the sigma meson 



tions here. In ref. 13 



mass and the pion mass, respectively, so that he concludes that the NG theorem is violated 
at finite temperature. As discussed before, we do not regard these values as the physical 
meson masses and rather consider them simply as variational parameters. We calculate 
the sigma meson mass according to eq. (|75D instead. 

When we evaluate the sigma meson mass in the chiral limit in this approximation, how- 
ever, we encounter a difficulty due to the infrared singularity. As shown in Subsec. 



there is a factor P2['^] in the equation of the sigma meson mass eq. ([75|) . Its temperature- 
dependent part P2.t[^] has the form 



47r2 



dx- 



X 



X' 



1 + 



t) J e 



-(f) _1 



(79) 



We naively suppose that this temperature dependent part vanishes at T — > 0. This is 
indeed so for m ^ 0. However, this factor gives a non-zero contribution if m approaches 
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to zero for T — > keeping m/T finite. This phenomenon actually occurs in the strict 
chiral limit if we neglect quantum corrections, because one of the solutions of the SD- 
equation, Mx-, vanishes at T = 0. The fact that Mx/T approaches to a non-zero value at 
T = can be confirmed by numerical calculation. Thus the sigma meson mass m„{T — > 0) 
deviates from m„{T = 0) = 600 MeV. We consider that this is because we have taken only 
the temperature dependent part in P2['^]- In fact, in the chiral limit, the temperature- 
independent part P2,q[^] is also infrared divergent at T = and is to cancel the infrared 
singularity in P2,t["^]- So it may not be possible to separate the temperature-dependent 
part and the temperature- independent one in P2['^]- 

There are two possibilities to solve this problem. One is to introduce the explicit 
chiral symmetry breaking term —ca, as seen above. This makes Mx finite at T = 0, 
corresponding to the finite-mass pions. So if we take the pion mass to be considerably large, 
e.g., m^(r = 0) = 138 MeV, this problem apparently disappears. Another possibility, 
which we think is better, is to include quantum corrections, which makes Mx finite even 
in the chiral limit. The study along this line is given in the following subsection. 

5.2.2 Thermal and Quantum Corrections 

In this subsection, we include the quantum corrections following the renormalization tech- 
nique of the auxiliary field discussed in Sec. |^. We have then another parameter /i which is 
introduced in the dimensional regularization and corresponds to the renormalization scale. 
As /i is a free parameter, we must choose a suitable value /i at which the spontaneous 
symmetry breaking actually occurs. All the other parameters are determined at the cho- 
sen n as reasonable values. Before proceeding to finite temperature, we perform a simple 
estimation in our model. 

In QCD, the following relation is well satisfied. 

f^ml = 2m,{qq) (80) 

where q denotes the light u, d quarks. This relation is called the Gell-Mann-Oakes- 
Renner(GMOR) relation and is derived from the low-energy theorem of QCD. In the chiral 
perturbation theory which is one of the low-energy effective models of QCD, this relation 
is also satisfied to the first order in the quark mass. This equation shows that the pion 
mass square is proportional to the quark mass. The deviation from it corresponds to higher 
order effects of the chiral perturbation theory. Since the chiral perturbation theory is a 
perturbation expansion in both external momenta of the NG bosons and the current quark 
masses, if this relation is satisfied up to the relatively large quark mass, it means that the 
chiral perturbation theory is valid up to this momentum scale. We recently investigated 
the behavior of the GMOR relation by constructing the pion as a relativistic bound state 
of quark and antiquark pair using the Bethe-Salpeter equation ||22|. We showed that the 
GMOR relation is quite well satisfied at least up to the strange quark mass (~ 150 MeV) 
region. 

In our linear sigma model the relation corresponding to GMOR relation is given by 
eq. ( |7BD . We can evaluate the c-dependence of in the following way. First we determine 
the values of c and the other parameters at m,r = 138 MeV. Then by changing the value 
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of c with the other parameters fixed, we calculate Mxy and Mx, and the meson masses. 
The obtained results are shown in Fig. |^. We have confirmed that these results are 
almost independent of the renormalization scale. These three values show almost linear 
dependence on c. In the chiral limit (c = 0), the sigma condensate which corresponds to 
the pion decay constant is about 87 MeV. This is very similar to the one obtained in the 
chiral perturbation theory p3[| which gives 88 MeV. The sigma meson is about 490 MeV 
(we have fixed rria = 600 MeV at m^r = 138 MeV) and the pion mass of course vanishes in 
the chiral limit. Thus we see that the linear sigma model gives consistent results with the 
chiral perturbation theory. 

Now we proceed to finite temperature. Once the quantum corrections are included, 
Mxy and Mx do not agree with the sigma and pion masses even at T = 0. We need to 
determine them by solving the SD-equations and the gap equation for the sigma condensate. 
We show the results of the initial parameters obtained in this way in Table. |l|-0. The value 

is the point where the uo curve crosses the temperature axis. We see that for smaller 
/i, Mx can be negative so that the effective potential becomes complex unphysical. For 
larger becomes positive so that the spontaneous symmetry breaking does not occur. 

As a result a proper range of fi is naturally determined as is given in the table. Because 
the parameter fi is free within this range, we show as an example a result for = 320 MeV 
in the case of m^{T = 0) = 600 MeV. 

We show the temperature dependence of the solutions of the SD-equations, the physical 
meson masses and the effective potential in the chiral limit in Figs. Temperature 
dependence of the sigma condensate is also shown in Fig. ^. The effective potential shows 
the signature of the first order phase transition, whose behavior is similar to the case 
without quantum corrections. 

We see also from the solutions of the SD-equations and the sigma condensate that the 
equations have two solutions at some temperature. This is a typical feature of the first 
order phase transition, that is, the upper solution corresponds to the (local) minimum of 
the potential and the lower to the maximum of the potential. 

As the temperature increases from zero, the sigma condensate decreases, jumps to zero 
at Tc and remains zero above T^. Likewise, in the figures Mxy [Mx) changes along the 
solid (dashed) line, jumps to the lower solid line at Tc and increases along this line. 

If the deconfinement transition takes place at the same critical temperature as the 
chiral phase transition, there are free quarks and gluons for T > Tc. But recent lattice 
simulations suggest that there are some correlations of quarks and gluons above p4 |. 
An instanton liquid model also favors this result This may imply there exist hadronic 
excitations even above Tf.. 

The fact that the phase transition is the first order agrees with other mean field 
approaches]^, |§, 0- It is, however, generally believed mainly from the analysis of the 
renormalization group equation that the actual phase transition of the 0(4) linear sigma 
model should be of the second order. We comment on this problem in the following section. 

The critical temperature is about 184 MeV for fi = 320 MeV and about 178 MeV for 
/i = 500 MeV. For the case of m^(T = 0) = 1 GeV, it is about 220 MeV for fi = 400 MeV 
and about 230 MeV for /i = 580 MeV. So the renormalization scale affects the critical 
temperature little. These values are reasonable since other approaches in the linear sigma 
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model predict similar values. 

We note that there is no peculiar behavior in the curves of the sigma meson mass, 
shown in Fig. ^. This is because Mx does not vanish even at T = due to the quantum 
corrections, shown in Fig. |^. On the other hand, the pions remain massless below the 
critical temperature even in this case. 

Finally we show the cases with the explicitly chiral symmetry breaking, i.e. m^(T = 
0) = 138 MeV in Figs. ^ 



6 Comments on Effects of Other Loops 

We have so far calculated the physical quantities in the Hartree-Fock approximation. This 
means that the solutions of the SD-equations for the propagator, Mxy and Mx, are 
momentum-independent. In the super-daisy diagrams, this assumption is quite natural 
because no external momentum enters into loops. In fact many authors use this approxi- 
mation within the formalism of the super-daisy diagrams in the mean field approach]^ ^ 0. 

However, note that near the critical temperature there is no reason to believe that 
only the super-daisy diagrams are dominant. Within the framework of the Hartree-Fock 
approximation, the first-order phase transition takes place as in our case. (If we take the 
large limit in the CJT approach, the result is the second-order, see ref.[|13|). But the 
analyses of the renormalization group equations etc. suggest that the phase transition of 
the 0(4) linear sigma model is the second-order |]26|, Arnold and Espinosa pointed out 



that other loop diagrams than the super-daisy diagrams are important near the critical 
temperature So the result for the order of the transition cannot be trusted in the 

calculation with the super-daisy diagrams only. 

The next leading loop diagram contributing to the effective potential is a setting- 
sun(SS) type diagram. This diagram contributes to the two-point Green function as in 
Fig. |T0|. 

It is obvious that the external momentum enters into the loop and therefore the self- 
energy depends on the external momentum. In our CJT formalism, this means that Mxy 
and Mx become momentum- dependent. Furthermore, we may have to take into account 
the wave function renormalization by this diagram. As far as we know, there is no cal- 
culation in which this diagram is incorporated consistently. Here we briefly estimate the 
contribution of the SS diagram with an approximation in order to see whether the order 
of the transition changes or not. The term of the effective potential corresponding to the 
SS diagram is given by 



4A^ 



Wei i {2GaMGhd{q)GU-p-q)+Gac{p)GM{q)GM{-p-q)} (81) 



Hereafter we consider a simple case with A^ = 1 since we are now interested in the phase 
transition rather than the hadron phenomenology. The CJT effective potential is now 
given by 

m,^) = lm'<Pl + ^<Pt + U \nG-\k) + l4 V-\k)G{k) 
2 24 2Tp,k 2Tf3,k 
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+^(^^^ G{k))' - ^<Plf^ ^ ^ G{p)G{q)G{-p - q) (82) 



with 



V-\k) = e + m^ + ^(g (83) 



From the extremal condition, we obtain the SD-equation for G. 

G-\p)=V-\p) + \i G{k)-^<Pli G{k)G{-p-k) (84) 
2Tp,k 2 Ti3,k 

Since the integral of the last term in eq.(0) depends on the external momentum, it is 
natural to put the form of G as 

G{p) = ^ ^ff^ (85) 

Here we take a very simple form in which any external momentum dependence in the 
parameters are neglected, i.e., A = \ and M(p) = M. It can be calculated relatively easily 
in the CJT formalism by this procedure. Though this is rather a strong approximation and 
seems to be similar to the super-daisy approximation, we see that the order of the phase 
transition can change in some case. 

We have calculated two cases which neglect quantum corrections. One corresponds 
to the small A, A ~ 0(1), and the other the large one, A ~ O(IOO). The latter is the 
order used in the low-energy effective model of the vr, a mesons. The results are as follows. 
Without the SS diagram, if the coupling constant A is small, the weak first order phase 
transition occurs. We see that in this case the order becomes second by including the SS 
diagram, while in the strong coupling case the order does not change. So the contribution 
of the SS diagram in this approximation changes the order of the phase transition from 
the first to the second if the original first order is rather weak. 

Chiku and Hatsuda calculated the SS diagram in the optimized perturbation theory [|T0| . 
They include the resummed mass parameter in the lagrangian from the beginning and 
renormalize it first with the approximation that the mass parameter is momentum inde- 
pendent constant factor. We think that their approximation is at the same level as ours 
used above. Nevertheless they suggested that the SS diagram may change the order of the 
phase transition from the first to the second |TT[]. 



Nachbagauer also calculated the SS diagram in the one-component 0^ theory in the 
CJT action]^. He calculated the self-energy in a certain approximation in order to get 
rid of the UV-divergence and evaluated the Green function in the Fade approximation. He 
also pointed out that this diagram is important near the critical temperature, although he 
does not evaluate the order of the transition. 

In conclusion, it seems to be necessary to include the effects of the higher-order dia- 
grams, if we discuss the behavior near the critical temperature. Further investigation is 
needed in the future. 
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7 Summary and Conclusion 



We have studied the 0{N) symmetric hnear sigma model at finite temperature in the 
Hartree-Fock approximation using the CJT effective potentiaL This model has been con- 
sidered as the low-energy effective model of the low-energy mesons, i.e., the sigma meson 
and the pions. It has so far been understood that the NG theorem is not satisfied, or the 
pions are massive at finite temperature in the CJT formalism unless the large limit is 
taken. We showed that this is caused by an improper derivation of the effective potential, 
such as an 0{N) violating form for the solution of the SD-equation for the propagator. If 
0{N) symmetric form is chosen and one defines the meson masses as the second derivatives 
of the effective potential, the NG theorem is always satisfied even at finite values of A^. 
This is also satisfied at any order of the loop expansion of which is the contribution of 
the two-particle irreducible vacuum graphs. 

We have also used a renormalization prescription of the CJT effective potential from the 
analogy to renormalization of the auxiliary field. This renormalization works apparently 
both in the symmetry broken and restored phases while the conventional renormalization 
works only in the symmetry restored phase. 

When one is interested in effects of the thermal contributions at finite temperature, 
quantum corrections are often neglected. The classical equation of motions are satisfied at 
T = in this case and the thermal effects are added as temperature increases. We have 
solved the gap equation for the sigma condensate, the SD-equations and the meson masses 
numerically in this approximation. As a result the first order phase transition occurs, 
because this corresponds to the mean field approach. Related to this approximation, we 
encounter a problem in calculating the sigma meson mass in the chiral limit because of 
the infrared singularity. This difficulty is removed if one includes the quantum corrections 
because the solutions of the SD-equations which cause the infrared singularity remain finite 
in the chiral limit at T = 0. So consistent treatments including the quantum corrections 



which is discussed in Sec. ^ are needed to obtain more realistic results. In Subsec. p. 2. 2 
based on the above notion, we incorporated the quantum corrections. As a result the 
infrared singularity disappeared and the behavior apparently became better for both the 
chiral limit and the finite pion mass. 

The resulting phase transition in 0(4) linear sigma model is of the first order which is 
consistent with other mean field approaches. However the second order phase transition 
is reported in this model in the renormalization group equation analyses. In fact near the 
critical temperature other loop contributions may become significant and we showed the 
simple estimation of the effects of the setting-sun diagram in the CJT approach. As a 
result the order of the phase transition can change with this inclusion. Therefore further 
investigation is needed if one studies the order of the phase transition more rigorously. 

Appendix 
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A The 0{N) linear sigma model with auxiliary fields 



We applied the technique of the renormahzation of auxihary fields for the CJT approach 
in the main text. In this appendix we analyze the 0{N) linear sigma model at finite 
temperature by using auxiliary fields in order to show a similarity with the CJT approach. 
We also point out that if we analyze the 0{N) linear sigma model at finite temperature with 
the auxiliary fields from the beginning, some difficulties associated with the NG theorem 
arise again. Such an approach was carried out in ref. pO[, where one auxiliary field was 



introduced. However this causes the violation of the NG theorem unless some 1-loop terms 
are added in the sigma and pion self-energy. We require a more general form of auxiliary 
fields which do not break the 0{N) symmetry in order to relate it to the CJT approach. 

For the original lagrangian, after shifting the fields 0" + 0°, we introduce an 

auxiliary field with an N x N symmetric matrix. 



exp 



j Vx exp l~ J d'^^ 



3N 



X 



2Ac 
?,N' 



3N 



(86) 



We take the saddle point approximation for x""^, i-e-; x"^ ~ Xc^) ^-^id then the effective 
potential at the 1-loop order is obtained. 



2 , ^ r^2^2 



3A^ 



8A(A^-f8) 



Xc 



SN^"" N + 8 3N^^^^ N + 8 



X 



ab 



^7) 



From the extremal conditions for Xc 'Pt^ obtain the following equations. 



= + — + 



2A 2A 1 



+ 3 ^ {N + A)/N \ 



3iV ?> N + S%,k\k'^ + Ml k^ + Ml j 



Ml 



2 2A 2A 1 



(iV + 4)(iV- l)/(iV + 8) 



1 



Tp,k + Ml 



+ 



where we defined a and /5 by X, 

Ml 



2A , 



A;2 + Ml 
2X{N -I] 



+ 



A;2 + Ml 



1 







ab , a<Pi<f>l 



(89) 
(90) 



3A^ 



+ 



3A^ ip,k fc2 + Ml J 
without loss of generality of x"^ and put 

{N + 2)a + p 



N + 8 



(91) 
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r2 2 , 2A . , (iV + 2)a + 3/? 



= m' + —d)i + ^ (92) 

The equations (|88|), (|8^) and (l90| ) are similar to eqs. (|3^) , (0) and (|3^) in the CJT approach, 



respectively. In fact the renormalized form of the effective potential ( pTl) can be obtained 
by requiring renormalization conditions for Mi and M2 at zero temperature, which means 
the renormalization of the auxiliary fields. Our renormaliztion prescription in the CJT 
approach is based on this idea. 

A difficulty occurs, however, when we consider the case of finite temperature. Our 
numerical calculation shows that the equations (|88|), (p9|) and (pO| ) do not have any physical 
solutions at finite temperature in contrast to the CJT approach. This means that the 
summation of the super-daisy diagrams does not work well and we may require some 
additional terms to recover it as in ref. [^. The CJT approach, on the other hand, does 
not have such a problem, though the forms of the equations are similar. We think that the 
CJT case corresponds to special auxiliary fields other than those of this section, but the 
relations between them remain as the future problem. 
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Table 1: Initial Parameters for T = 

m^ = MeV 



[MeV] 


fi [MeV] 


A 


Mx [MeV] 


MxY [MeV] 


[MeV^] 


[MeV] 


600 


320 


90.2 


41.5 


510 


-122375 


128 


600 


400 


102 


119 


542 


-94208 


105 


600 


500 


122 


183 


594 


-47679 


68.4 


600 


567 


139 


223 


632 








600 


TREE 


125 





600 


-180000 


132 


1000 


410 


156 


29 


670 


-218034 


130 


1000 


500 


177 


156 


715 


-128224 


93.2 


1000 


599 


208 


239 


774 








1000 


TREE 


347 





1000 


-500000 


132 



Table 2: Initial Parameters for T = (various values of /i) 



= 138 MeV 



[MeV] 


/i [MeV] 


A 


Mx [MeV] 


MxY [MeV] 


[MeV] 


600 


220 


66.9 


55.9 


460 


-96324 


600 


300 


77.8 


114 


493 


-83362 


600 


400 


94.7 


165 


540 


-59939 


600 


520 


119 


223 


603 


-5633 


600 


TREE 


118 


138 


600 


-151434 


1000 


380 


142 


54.8 


655 


-197043 


1000 


400 


146 


87.5 


664 


-180329 


1000 


500 


172 


184 


718 


-95776 


1000 


580 


196 


245 


765 


-106 


1000 


TREE 


340 


138 


1000 


-471434 



Table 3: Initial Parameters for T — (various pion masses) 
= 600 MeV, /i = 320 MeV 



[MeV] 


A 


Mx [McV] 


MxY [MeV] 


[McV] 





90.2 


41.5 


510 


-122375 


20 


89.8 


44.7 


509 


-121078 


60 


87.8 


64.8 


507 


-112237 


100 


84.6 


93.5 


504 


-97658 


138 


80.9 


125 


502 


-79760 
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Figure 1: Two-particle irreducible diagrams contributing to r2[4>c, G] up to the three-loop 
level. 




Figure 2: Two-particle reducible diagrams which do not contribute to r2[0c, G] up to the 
three-loop level. 
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Figure 3: The solutions of SD-eqs. for rricriT 
solid line is Mxy and the dashed one is Mx 
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Figure 4: for m^{T = 0) = 600 MeV and // = 320 MeV 
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Figure 5: Effective potential for mfj{T = 0) = 600 MeV and /i = 320 MeV. Solid lines 
denote the case of T =198 MeV, 194 MeV, ■ ■ -,174 MeV from top to bottom. 
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Figure 6: cTq, and for /i = 320 MeV 
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Figure 7: The SD-solutions m<,(T = 0) = 600 MeV, m^{T = 0) 
MeV. The sohd hne is Mxy and the dashed one is Mx 



138 MeV and ii = 320 
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Figure 9: Sigma condensates (Tq for 'm„{T = 0) = 600 MeV, n 
0) = 0, 20, 60, 100, 138 MeV from left to right. 



320 MeV and m^{T 




Figure 10: The loop diagram contributing to the two-point Green function which comes 
from the setting-sun type diagram. 
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